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Abstract 

The collective oscillations of ID repulsive Bose gas with external harmonic confinement in two 
different regimes are studied. The first regime is the mean field regime when the density is high. The 
second regime is the Tonks- Girardeau regime when the density is low. We investigate the resonances 
under periodic modulations of the trap potential and the effective nonlinearity. Modulations of 
the effective nonlinear coefficient result from modulations of the atomic scattering length by the 
Feshbach resonance method or variations of the transverse trap frequency. In the mean field regime 
we predict the bistability in the nonlinear oscillations of the condensate. In the Tonks- Girardeau 
regime the resonance has the character of a linear parametric resonance. In the case of rapid strong 
modulations of the nonlinear coefficient we find analytical expressions for the nonlinearity managed 
soliton width and the frequency of the slow secondary oscillations near the fixed point. We confirm 
the analytical predictions by direct numerical simulations of ID Gross-Pitaevskii equation and the 
effective nonlinear Schrodinger equation with quintic nonlinearity and trap potential. 

PACS numbers: 03.75.Kk, 42.65.-k, 42.50.Ar 
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I. INTRODUCTION 



Low dimensional Bose-Einstein condensates (BECs) in highly asymmetric traps have 
recently been achieved, which opens new possibilities in the investigation of clouds of bosonic 
and fermionic gases Low- dimensional bosonic systems have many remarkable 

properties which distinguish them from 3D systems. One of them is the growth of the 
interaction when the density is decreased. As a result the system enters into the Tonks- 
Girardeau (TG) regime. The properties of the Bose gas then coincide with the gas of free 
fermions. This follows from the exact solution to theproblem of hard-core bosons with 
repulsive interaction obtained by Lieb and Liniger |5[ . The opposite regime with weak 
interaction (high density case) is the mean field (MF) regime Qj. 

The TG and MF regimes can be characterized by the parameter 7 which is equal to the 
ratio of the interaction energy and the kinetic energy of the ground state of gas, i.e. 7 = 
^giD / (f^'^iT'iD) ■ Here m is the atomic mass, Qid is the one-dimensional coupling constant, and 
riiD is the ID density. The case 7 <^ 1 corresponds to high densities, when the description 
by the mean field theory is valid. The case 7^1 corresponds to the strong repulsive 
interaction - the Tonks- Girardeau regime. The modern experiments with Bose gas in highly 
elongated traps now are in the region 7 ~ 1 

One of the important phenomena for the experiments is the collective oscillations of the 
Bose gas in the different regimes. It is particularly interesting to investigate the dynamics 
of breathing and dipole modes for MF and TG regimes, as well as the crossover between 
them. The theoretical predictions for the frequencies in harmonic longitudinal traps V{z) = 



mu;^2;^/2 are obtained in Refs. \8L lOlllO 
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12| . It is shown that the frequency of oscillations 



in the mean field regime is y/Suz, while it is 2ujz in the TG regime. The last value coincides 
with the one observed for the thermal gas. For the 3D cigar the frequency is ^/hjluj^. 
These investigations involve the combination of the exact solution by Lieb and Liniger, 
the local density approximation, the hydrodynamic equations, and the extended nonlinear 
Schrodinger equation. Note that information for oscillations in the high-dimensional regimes 



can be found in the review 



13|. In this work we shall study resonances in the oscillations of 



ID Bose gas in both regimes. The periodic and random modulations of the trap potential 
and the atomic scattering length are subject to our analysis. 

The MF regime is described by the ID Gross-Pitaevskii (GP) equation with two-body 
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interaction 

iHt = -^^zz + V{z, t)ij + giDT{t) I V-l V, (1) 

This equation is derived from the 3D GP equation in a strongly anisotropic external poten- 
tial. The dynamics in the radial direction is then averaged out 0] and the longitudinal profile 
of the wave function satisfies Eq. (P). The wave function is normalized so that the number 
of atoms in the BEC is N = J \tp{t,x)\'^dx. V{z,t) is the longitudinal trapping potential, 
which is assumed in this work to be harmonic V{z,t) = mujlz^ F {t) / 2 . The function F{t) 
describes the variation in time of the trap. The effective nonlinear coupling constant is gio- 
In the case of a harmonic transverse trap potential muj\{x'^ + y^)/'2,, Eq. (Q) is valid under 
the assumption uj±_ ^ uj^ and we have qid = 2hasUJ±. The function r{t) describes the vari- 
ation in time of the effective nonlinearity. A first method to vary the effective nonlinearity 
is to modulate in time the transverse trap width or equivalently the transverse frequency 
The modulations of the transverse frequency impose variations of the BEC density in 
the (x, yl-plane, which in turn involve variations of the nonlinear interaction il4|l . Nonlinear 

nn 

resonances in 2D BEC for such modulations have been studied in |15l.llfil|. A second method 
to vary the effective nonlinearity is to modulate the atomic scattering length by the so-called 
Feshbach resonance technique [l7|. 

Theoretical and experimental studies have demonstrated that variation of the s-wave 
scattering length, including a possibility to change its sign, can be achieved by using the 
Feshbach resonance 

Qs is the value of the scattering length far from resonance, B{t) is the time-dependent 
external magnetic field, A is the width of the resonance, and Bq is the resonant value of 
the magnetic field. Feshbach resonances have been observed in ^^Na at 853 and 907 G jl7|. 
in ^Li at 725 G and in ^^Rb at 164 G with A = 11 G 0. In the case of resonance 
dynamics where is slowly varying and keeps a constant sign atom losses are negligible. 
However atom losses may be important when crossing the resonance jl7l.l2fl|. This is the case 
of the ^^Na condensate where it is necessary to cross the Feshbach resonance to change the 
sign of the atomic scattering length. The approach developed in this paper should then be 
modified to take into account this phenomenon. If we are not close to the resonance, losses 
are small and in a first approximation they can be taken into account by a time-varying 
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number of atoms in the effective variational equation for the width. This should lead to 
damped secondary oscillations with a time scale larger than the oscillations time. However 
atoms losses can be minimized down to a negligible level by certain experimental control 



that have been implemented in particular in the *^Rb case 



Ref. 
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2J]. Furthermore in 



18[ it is demonstrated in the case of ^Li that a change of the sign of the scattering length 



can be obtained without crossing the resonance by the so-called coupled channel method. 
Our study of the nonlinear management is triggered by these experimental achievements. 

According to Kolomeisky et al. [2^ the TG regime is described by the nonlinear 
Schrodinger equation with quintic nonlinearity and trap potential 

iHt = -—(P,, + V{z,t)(t) + ^\(t>\% (2) 

It is known that this model does not capture every aspect of the dynamics of an atomic gas 
in the TG regime, in particular it overestimates the coherence in interference patterns at 
small number of particles [2^. However Eq. Q has been shown to reproduce the collective 



spectrum of a gas in the TG regime within a local density approximation [11|, |27|. Fur- 
thermore we shall show that it also gives the correct frequency of oscillations as predicted 
theoretically and observed experimentally 0]. As our attention is focused on resonance phe- 
nomena, Eq. Q seems to be a good model, both physically relevant and mathematically 
tractable. Using the time-dependent variational method or the hydrodynamic approach we 
shall derive the equation for the TG gas width and study resonances in oscillations under 
periodic and random modulations of the trap potential. 

The paper is organized as follows. In Sections II-IV we analyze the nonlinear resonances 
in gas oscillations in the mean field regime using the ID mean-field GP equation. We apply 
a time-dependent variational approach and introduce action-angle variables in Section II. 
Section HI, resp. IV, is devoted to the resonances driven by periodic and random mod- 
ulations of the trap frequency, resp. the nonlinearity. In Section V we address the same 
problems for the ID TG regime. 
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II. THE MEAN FIELD CASE 



A. Variational approach 

We first put Eq. ((T)) into dimensionless form by setting t' = u^t, x = z/l^, h = a/V(™^) 
and u = ^y2\as\uj±/oOz'^- In tlie following we omit primes, so the mean-field GP equation 
reads 

1 1 ~ 

iut + -Uxx - -F{t)x^u - V{t)\u\^u = (3) 

where F(t) = Fit/uj^) and r(t) = r(t/ci;^)sgn(as). We apply a variational approach using 
the Gaussian ansatz 

«(t,x) = Aexp(^-^-^-z^^ icp{t)\ (4) 

Note that the equation for the gas center of mass is decoupled from the equations for 
oscillations so we did not take it into consideration. The ansatz yields a closed-form evolution 
equation for the width 

au = -,- f{t)a + ^ (5) 

where f{t) = F{t), -f{t) = Pt{t), P = J \u\'^dx/V2^ = {2\a,\u^N)/{V2^ujJ,). In the 
repulsive case > and in absence of modulations F = 1, F = 1, we have / = 1 and 
"f = P = N/N* with A^* = (^/27iuJzlz) /{'2asUJ±). If the number of atoms is large enough 
N ^ N* then P ^ 1 so that we can neglect in Eq. (jS)) and the fixed point is given by 

S = P'^' (6) 

In the dimensional variables we get Lc = [{V2aslluj±N){y/Truj;,)Y^^, that agrees (up to a 
numerical multiplicative constant) with the Thomas-Fermi value for the ID BEG width. 



B. Action-angle variables 

We assume in this section that 7(t) = 7o + 7(t) where the average nonlinear coefficient 
7o > which corresponds to the repulsive case, while 7 represents a zero-mean periodic or 
random component. Similarly, we take f{t) = 1 + f{t). The unperturbed problem consists 
in taking 7(t) = f{t) = 0. Assume that 70 is large so that the kinetic term can be 
neglected in Eq. The energy E of the unperturbed BEG is given by: 

m = ^a?(t) + U{a{t)), U{a) = + ^. (7) 
5 



In absence of fluctuations the energy E is an integral of motion. The BEC width obeys a 
simple dynamics with Hamiltonian structure 
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(8) 



with q = a, p = at- The potential U possesses a unique minimum Ug = 70^''^ which is a 
stable flxed point with oscillation frequency u = ^/3. The corresponding ground state has 
energy Eg = U{ag) = {?,/2)^T. 

If the initial conditions (a(0),af(0)) correspond to an energy above Eg, then the orbit of 
the motion is closed, corresponding to periodic oscillations. In order to explicit the periodic 
structure of the variables a and at-, we introduce the action-angle variables. The orbits are 
determined by the energy imposed by the initial conditions 

E = ia2(0) + f/(a(0)). 

For E > Eg, we introduce ai{E) < a2{E) the extremities of the orbit of a for the energy 
E. They are the positive solutions of the cubic equation U{a) = E and they are given by 
(J = l,2) 

1/2 



aj{E) 



2E 



arccos 



E 



cos 

3/2' 



-lF27r 



The action I is deflned as a function of the energy E by 

If 1 r^-iiE) 



T{E) = — fpdq 



^/2E - 2U{h)dh. 



(9) 



ai(E) 



The motion described by (jH)) is periodic, with period 

J P Ja,{E) ^2E - 2U{b) 

or else T{E) = 2ti^{E). The angle is deflned as a function of E and a by 



(10) 



HE, a) 



" dp , 2tt 
m^' = -T{E) 



dh 



^2E - 2U{h) 

The transformation {E,a) {1,4') can be inverted to give the functions and A{I,(f)). 
The BEC width oscillates between the minimum value ai{E) and the maximum value a2{E). 



The energy E as well as the action / are constant and determined by the initial conditions, 
so the evolution of the BEC width is governed by 



a{t)=AmE),4>{t)), 

m = 0(0) - 

For E close to the ground state energy Eg, we have 



a,(E) = a, + (-l)y^(E-E,), J = 1,2, (11) 

T(^) = ^, I{E) = ^~^\ AiI,<P) = ag + V2Icos{^). (12) 
V 3 V 3 



For large energies E ^ Eg, we have 



2,E / 3E 

ai{E)=agj^, a2{E) = agJ — , (13) 

nE) = n, I{E) = ^, ^(/,0) = v^Vl + cos(0). (14) 



III. RESONANCES IN THE MEAN FIELD CASE DRIVEN BY A TIME- 
VARYING NONLINEARITY 

In this section we address the role of a time- varying nonlinearity and assume that the trap 
is stationary f{t) = 0. We shall mainly focus our attention to the periodic management 
r(t) = Fq + Fisin(r2gt), but we shall also consider random fluctuations of the effective 
nonlinearity in the ID GP equation. 



A. High-frequency nonlinear management 

We shall first address the case where the oscillation frequency of the nonlinear manage- 
ment is much higher than the trapping frequency, i.e. Qg ^ Uz- We must also assume 

^ to prevent from exciting the transverse modes. In such a case the influence of 
the nonlinear management is negligible unless the nonlinear management amplitude is large. 



The problem of nonlinearity management for BEC, the so-called 
management, has already been considered for ID BEC in Refs. 




in Refs. j29|, |30|. In Ref. [2^ the authors were the first ones to propose a technique of FR 



eshbach resonance (FR) 



3,Q 



and for 2D BEC 
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management, based on a time-periodic change of the magnitude and sign of the scattering 
length by a resonantly tuned ac magnetic field. The FR management resembles the so-called 
dispersion-management (DM) technique in fiber optics, which is based on a periodic con- 
catenation of fibers with opposite signs of the group-velocity dispersion. The DM technique 
has been shown to support robust breathing pulses, the so-called DM solitons celebrated in 
optics jsil. The FR technique is shown in Refs. [24, 28 1 to drive stable localized structures, 
named FR managed matter wave soliton. Here, in distinction from 2^ 2^, we give the an- 
alytical expressions for the fixed point (corresponding to the FR managed soliton) and the 
frequency of slow secondary oscillations. More precisely, we introduce the small parameter 
6 = uJz/Qg and we assume that the nonlinear management amplitude is large compared to 
the average value. We write accordingly 

lit) = 7o + y sin(^) 

and perform an asymptotic analysis 5 — > following the Kapitsa averaging theorem. We 
expand a(t) = ao(t) + 6ai(t,t/6) We substitute this form into Eq. ^ and get a 

compatibility condition which reads 



— - ao + -2 + -5 , 



a: 







a, 







a; 



(15) 







while the first-order corrective term can be expressed as 



ai(t,r) 



71 
alit) 



sm r 



Note that the initial conditions for the slowly varying envelope are 



ao(0) = a(0), a'o(O) = a'(0) - 



71 



a(0)2' 

where a(0) and a'(0) are the initial values of the width and its time-derivative. The ground 
state can then be analytically studied. Let us first deal with the case 70 = 0. We introduce 



the critical value 7c = ^4/27 ~ 0.62. It is found out that, if |7i| < 7c, then Eq. (|THjl admits 
a unique fixed point describing the width of the ground state 



cos 



-arccos 
3 



7? 
7c^ 



If |7i| > 7c, then 



a: = 
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FIG. 1: BEC width (picture a) and oscillation frequency (picture b) as predicted by the theoretical 

1/3 

model (in dimensionless units). For 70 ^ 1 the BEC width is 7q and the oscillation frequency is 

Qg increases with 71 and goes from the value 1 for 71 = to the asymptotic behavior 
ag ~ 'fi^^ for large 71. The linear stability analysis of the effective equation ()15|) shows that 
the fixed point is stable. If the initial condition is close to this point, then the width a 
oscillates around the value ag with the oscillation frequency u = — which increases 
from 2 for 71 = to its limit value \/6- In the general case 70 7^ 0, 71 7^ 0, there exists a 
unique fixed point which is the unique positive zero of the equation — + 790'^ + 7^ = 0. 
It is plotted in Figure [T^, and the corresponding oscillation frequency is plotted in Figure 

Hb. 

We have carried out numerical simulations of the ODE model (0) and the PDE model Q 
to check the theoretical predictions of the asymptotic analysis 5-^0 (see Figure E]). Note 
that stable BECs can be achieved with a negative (i.e. attractive) or positive (i.e. repulsive) 
average nonlinear coefficient 70. 

Let us estimate the parameters for a realistic experiment. For the ®^Rb condensate with 
as = -0.5nm, u;^ = 27r • 360Hz, 00^ = 27t ■ 14.4Hz, = 10^ I, = 100 /zm, Qg = lOu,, Fq = 1, 
Fi = 2, we find that the condensate width is ag ~ 0.71^ and the frequency of the secondary 
oscillations is a; ~ 2.6lJz. 
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FIG. 2: BEG width in presence of high-frequency periodic modulation of the nonUnear coefficient. 
Here 7(t) = 70 + 2sin(10t) with 70 = 1 (picture a) and 70 = —1 (picture b). The sohd and 
dashed hues stand for full numerical simulations of the PDE system with two different initial 
conditions (dashed lines: theoretical ground states), while the dotted lines represent the theoretical 
predictions. 

B. Resonant nonlinear management 

In this section we address the case where the periodic nonlinear management is resonant 
or close to resonant. We shall focus our attention to the particularly interesting case where 
the number of atoms is large, which in turn implies that the dimensionless parameter 70 
is large (say at least 5). We shall see that a periodic modulation of the nonlinear coefficient 

7(^) = 7o + 7i sm{nt) 

may dramatically modify the dynamics, and this phenomenon will be noticeable when the 
dimensionless parameter 71 is of order 1. If 70 > 5, then Eq. can be simplified into 

au = -a + ^. (16) 

In absence of periodic modulation 71 = 0, the ground state is Qg = 1^/70 and the oscillation 
frequency of the BEC is = y/S in dimensionless units. We now consider a periodic 
modulation of the nonlinear coefficient with an amplitude 71 smaller than 70 and with 
frequency fl of the same order as ur. We expand a = Ug + d where d satisfies 

d + uP'd + bo} + Aa^ = e sin(f2t) + sin(f2t)a, 
10 



]^ /3 2 /3 2 /3 1 

with 6 = — 37q , a = 47q , e = — 7i7o , and €2 = —2717,^ . This is the equation for a 
nonhnear oscillator with external and parametric drives. The analysis of this problem can 
be carried out by applying the standard method described in Ref. js^. We consider the 
harmonic expansion of a 

a = di cos(z/) + flo + cos(2z/), 

where z/(t) = Qt + 6{t). We substitute this ansatz into Eq. (fTHIl and we collect the terms 
with the same harmonic. This yields a system of differential equations that give the zeroth 
and second harmonics in terms of the first one, and a compatibility condition which reads 
as a system of two differential equations for the first-harmonic Si and the slow phase 9 

e2COS^(6')\ €28 (ld\ ecos(^^)ai\ Te^ai 



dtai 



sin(0),(17) 



2 - fi2 ^g^g g / e2pos2(0)\ cos(0) 



+ (-9^ + 8^ ) + (s?^ - 12^ J 

The first harmonic is dominant so that the BEC width evolution can be described at first 
order as a{t) = ag + di{t)cos{Qt + 0{t)). The stability analysis shows very interesting 
features. It appears that the strongest resonance occurs when the nonlinear management 
has the frequency 




n^ = uj\li + -[j^J , (19) 

which is above the oscillation frequency of the BEC. Resonance is still noticeable for a 
modulation frequency f2 in a vicinity of fic with bandwidth of order (71/70)^''^. More 
precisely, if < Qc, then there exists a unique fixed point to system ()17H18|) . The value of 
di corresponding to this fixed point is 



_ 1 

(^osc TT 



2771 + 3a/2472(u;2 _ ^2)3 + 81^2 



27o'/'(c^^-^^^ 



2771 + 3v/24^^(IJ^^TF)3:^81^ 



1/3- 



A linear stability analysis shows that this fixed point is stable. HQ > Qc, then there are three 
fixed points. A bistability in the condensate oscillations occurs. The intermediate fixed point 
is always unstable. The upper fixed point is the continuation of the fixed point exhibited 
in the case Q < Qc- It can be observed by increasing slowly and carefully the modulation 
frequency from a frequency below Qc to a frequency above Qc- The lower fixed point is the 
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FIG. 3: Oscillation amplitude as predicted by the theoretical model (solid lines). Here 70 = 10 
and 71 = 1. Comparisons with simulations with the ODE model (crosses) with the PDE model 
(circles). 

one that is observed when imposing without any particular precaution a modulation with 
frequency Q > Q^i or when perturbing the metastable upper fixed point. The value of the 
stable fixed point is 



osc 



2\/^7o^'(^^^-^^)cos(^ 



Cose = arccos 



7i| 

3/2 



Thus the system encounters a jump in the oscillation amplitude when the modulation fre- 
quency crosses the value flc- Evaluating the two expressions of aosc just above and just below 
the critical frequency Qc, "we get that the amplitude of the jump is aosd^c) ~ o^osd^t) — 
l7i|^/'. 

In Figure El we plot the values of the oscillation amplitude as a function of the frequency 
of the nonlinear management in the case 70 = 10 and 71 = 1. The theoretical prediction is 
that the modulation frequency driving the most resonant response is Qc — 1-823 which is 
above the eigenfrequency cu = -\/3 ^ 1.732. We have carried out numerical experiments with 
the ODE model (0) and the full PDE model Q to check this prediction. The results for the 
frequency Q = 1.5 are plotted in Figured For each frequency we can detect the amplitude of 
the oscillations and report on Figure El to compare with the theoretical predictions. We have 
numerically determined that the frequency driving the most important resonance is ~ 1.84 
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FIG. 4: BEC width versus time for a nonlinear management with frequency = 1.5. We compare 
the theoretical envelope of the oscillation with numerical simulations of the ODE model (picture 
a) and the PDE model (picture b). Here 70 = 10 and 71 = 1. For the simulation of the PDF model 
the initial state is the Gaussian ansatz with ag = 2.15, which corresponds to the theoretical fixed 
point in absence of nonlinear management. The solid lines represent the theoretical slowly varying 
envelope a{t) = ag + 2aosc\ sin(a;osc^/2)|. 

which is indeed above the eigenfrequency V3 and very close to the theoretical prediction f2c. 
The numerical simulations confirm the jump in the oscillation amplitude when crossing the 
critical frequency. 

Let us estimate parameters for a realistic experiment. The magnetic trap can be taken 
with parameters u)± = 27r 10'^ Hz, uJz = 2tt 10 Hz, and the number of atoms of ^^Rb 

= 5 10^. For the external field B = 159 G, = 0.8 nm (repulsive gas). Then, for 
Q = 18.25 Hz, we should observe large oscillations with aosc ~ 2.1/^, while at = 15 Hz 
we should observe much smaller oscillations with aosc ~ 0.3/^. Practically, if we start from 
the fixed point ag given by (0), we should observe oscillations of the BEC width with the 
maximum ag + 2aosc- 

C. Random fluctuations of the nonlinear coefficient 

We assume here that the scattering length or the transverse frequency ujj_ is randomly 
varying inducing random fluctuations of the nonlinear coefficient of the GP equation 



lit) = 7o + 7i^W, 
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where r]{t) is a normalized random noise with standard deviation of order 1 and 71 > 
represents the amphtude of the random fluctuations. This model also presents interest for 
the study of nonlinear management schemes in spatial optical solitons if the widths of the 
nonlinear layers in arrays of waveguides are randomly distributed 3^. We use the angle- 
action formalism introduced in Section IllBI In presence of perturbations, the motion of 
a is not purely oscillatory, because the energy and the action are slowly varying in time. 
We adopt the action-angle formalism, because it allows us to separate the fast scale of the 
locally periodic motion and the slow scale of the evolution of the action. We assume that 
7i is smaller than 70 and introduce the dimensionless parameter S = 71/70- Thus, after 
rescaling t = the action-angle variables satisfy the differential equations 

where h(I, cj)) = — jo/AU, </>) and uj(I) = ^.^7r\^ smooth functions and h is periodic with 
respect to (p with period 2n. Applying a standard diffusion- approximation theorem 
establishes that, for small 6, {I{T))r>o has the statistical distribution of a diffusion Markov 
process characterized by the self-adjoint infinitesimal generator 

where 

-1 p27r POO 

A{I) = - / h^{I, (l))h^{I, (j) - uj{I)t) iviOHt)) dtdcj). 
Jq Jo 

and the brackets stand for a statistical averaging with respect to the distribution of the 

noise rj. This means in particular that the probability density function of / satisfies the 

Fokker-Planck equation drP = Cp, pij = 0,1) = 6{I — Jq), where Jq is the initial value of 

the action at time 0. As long as the energy remains close to the the one of the ground state, 

we can use the asymptotic expansions (|12j) to expand the effective diffusion coefficient A{I). 

In the original time scale we then establish that the growth of the action is 

2 



where lu = ^/S, Wt is a standard Brownian motion (i.e. a Gaussian process with zero-mean 
and standard deviation v^), and 

POO 

ac= {ri{0)ri{t)) cos{ujt)dt. (20) 
Jo 
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If the BEC is in the ground state at time 0, then the BEC width a oscillates with the 
eigenfrequency to = y/S between the values a_ and a+ which evolve slowly as 



a±{t) = 7o ± 



1/3 , V2al7i 



2/3 

^7o 



Wt\. 



As a result, taking into account the periodic modulations, the rms BEC widthisrms2(a)(t) = 
'^clit/{uj'^%^^) where the root mean square (rms) is defined by rms(a) = ((a — (a))^)^^"^ . 

In the white noise case ('7(0)?7(t)) = 2cr^5(t), there exists a direct way to compute the 
growth of the rms amplitude. It consists in writing a closed form system for the second-order 



moments of the BEC width and its derivative. We expand a 



a and denote b = at. By 

considering the column vector X(t) = {^{a^) , ( 56 ) , ( 6^ ) ) , we get by applying Ito's calculus 
that X satisfies the closed-form system 

dX 



dt 



MoX + MiX + K), 



(21) 



where 



Mn 



V 



-Ijj 




2 




-2u^ 



Ml 






,,2^,-2 



y 8a,7f 7o"' j 






O 2 -4/3 

\2ac7i7o / 



uj = y/S, ac = cr^, Vq is the source term and Mi is a stochastic resonance term. It is negligible 
during the first steps of the dynamics, and becomes important only when the rms amplitude 
becomes of the order of ag = 7^'^. However, in that case, the linearization procedure is not 
valid anymore. Integrating Eq. (j2H) by neglecting Mi we get 

sin(2ti;t) \ 



rms {a){t) 



«c7i 

4/3 2 

7o ^ 



t - 



2u 



■ ■ 



(22) 



In Figure El we compare the statistical predictions with a set of numerical simulations of 
the PDE model. We have taken 7o = 7i = 10. We have adopted a stepwise constant 
model for ri{t). rj is constant over elementary intervals with duration and takes random 
values uniformly distributed between —1 and 1. With = 0.03 we have Oc = 0.005. The 



initial BEC width is ag 
noticeable. 



7o 



1/3 



2.15. The diffusive growth in agreement with Eq. ()22|1 is 
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FIG. 5: Rms of the BEG width in presence of random fluctuations of the nonhnear coefficient. The 
theoretical result [Eq. H22() ] is plotted in dashed line. The results of numerical simulations of the 
PDE model are plotted in solid line and correspond to the averaging of 1000 different realizations 
of the random noise r]. 

IV. RESONANCES IN THE MEAN FIELD CASE DRIVEN BY A TIME- 
VARYING POTENTIAL TRAP 

We first focus our attention to the periodic management f{t) = l + /i sin(f2t), and second 
we address random fluctuations of the trap. 

A. High-frequency periodic modulation of the trap 

We shall first address the case where the oscillation frequency of the trap modulation is 
higher than the trapping frequency, i.e. ^2^0; = -\/3. In such a case the influence of the 
modulation is negligible unless its amplitude is large. We introduce the small parameter 
S = 1/Q and assume that the trap modulation amplitude is large, of order 6~^. We write 
accordingly 



and perform an asymptotic analysis 6^0 following the same line as in Section IIII Al We 
get that the slowly varying envelope of the BEC width obeys the effective equation 



/(t) = l + ^sin(-) 




2n^ 
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with the initial conditions 

ao(0) = a(0), a[,(0) = a'(0) - 

The high-frequency modulation thus involves a shift of the effective potential which in turn 
imphes a shift of the ground state (see also 3^). 

B. Resonances due slow periodic variations of the trap potential 

The study is similar to the nonlinear management, and qualitatively the same conclusion 
holds true, especially concerning the bistable diagram. We shall only point out the main 
differences. First, the strongest resonance occurs when the periodic modulation of the trap 
has the frequency 

= c.y^i + i|A|2/3, (23) 

which is above the eigenfrequency of the BEC. Note that Qc does not depend on the nonlinear 
coefficient 70 and is proportional to the eigenfrequency of the BEC. The oscillation amplitude 
of the BEC turns out to be also proportional to the BEC width. The resonant bandwidth 
is of order l/ip/^cj. More precisely, if < fic, then the oscillation amplitude of the BEC is 



a. 



osc 



ttg 3 



27|/i|+3./24(a;2_l]2)3 + 8i/2 



27|/i|+3v/24(^2_^2)3 + 8ij2 



while for f2 > fi^ 



2W-(r]2_^2)cos(^ 



Clg V 3 V 3 



^osc = arccos 



1/1 



Evaluating the two expressions of aosc around the critical frequency Qc, we get that the 
amplitude of the jump is aosd^") - ciosci^t) = l/iT^^To^^- 

C. Random fluctuations of the trap 



We consider in this section a random modulation of the trap 

/(t) = l + /ir/(t) 
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with < /i ^ 1 and 77 is a normalized random noise. We once again use the action-angle 
formalism. We carry out an asymptotic analysis similar as the one presented in Section 
1111 C\ with the small parameter fi and the function h given by 

h{I,<P) = ^A'{I,<P)- (24) 

We then get the statistical distribution of the slow evolution of the action in terms of a 
Brownian motion 

where u = a/3 and ac is given by (pUj). As a consequence the BEC width a oscillates with 
the frequency uj between the values a_ and a+ which evolve slowly as 

a±{t)=ag±^f,j'J'\Wt\. 

UJ 

This means that the BEC spreads out at the diffusive rate rms^(a)(t) = acfi%^^t/{2uj'^). 
The doubling of the width is observed after a time of order uj'^/{acfl)- In the case of an 
optical trap imposed by a laser field whose intensity is fluctuating the typical fluctuations 
level is of the order of Uc ~ 0.01. If the trap frequency cu^ = 27r • 300Hz, then we predict that 
the doubling of the width should be observed after a time of the order of a few seconds. 



V. RESONANCES IN THE TONKS-GIRARDEAU REGIME 

As it is shown in Ref. 



25| the variety of properties of the hard-core Bose gas with repulsive 
interaction in the dilute regime can be described by the nonlinear Schrodinger equation with 
quintic nonlinearity 

fc2 2 fc2 

iHt = -^0.. + V{z, t)<P + ^101 V, (25) 

where V{z, t) is the time dependent trap potential. In principle it can include an anharmonic 
part together with the harmonic component. Below we will restrict ourselves the harmonic 
case V{z,t) = muj1z'^F{t) /2. The wave function is normalized to the number of atoms 
J = A^. This equation takes correctly into account the dependence on the density of 

;he energy of the ground state of ID Bose gas and reproduces correctly the collective modes 
It was shown in Q by means of numerical simulations for small number of atoms, 
that the interference effects are overestimated by this equation. It should be noted that 



18 



the nonlinear coefficient does not depend on the scattering length, i.e. the details of the 
interaction. Accordingly this parameter cannot be managed. Introducing x = z/lz-, t' = tUz, 
u = \/T^lz(j)/2^/^, Iz = \/ hj (muj) we can write the equation in the dimensionless form 

iut = -^u^x + ^f{t)x'^u + (26) 

where we have dropped the primes and f(t) = F{t/uJz)- As can be expected the change 
in the exponent of the nonlinear term induces strong differences between the MF regime 
and the TG regime. The quintic nonlinearity is especially interesting in the ID case as it 
represents the critical nonlinearity for the NLS equation. Indeed, for a given dimension d, 
there exists a critical exponent a = 2/d for the nonlinear term Iwp'^u which separates two 
different regimes jsTj. These regimes are different in terms of global existence, blow-up, 
instability growth,... In the ID case, this exponent is a = 2 (quintic NLS), in the 2D case 
it is o" = 1 (cubic NLS). In the BEG framework, this leads for the repulsive gas to the 
parametric instability in a trapped ID quintic system similarly to the parametric instability 
observed in a trapped 2D cubic system. For the attractive gas this leads to the collapse in 
a ID quintic system with a number of atoms > A^^^ as well as the collapse in a 2D cubic 
system for N > N^. 



A. Variational approach 

We apply the variational approach with the Gaussian ansatz (0)). We find that the 
equation for the phase is decoupled from the equation for the atomic cloud width that reads 

att + f{t)a = (27) 

where C = j \u\^dx = {'r/V2)N and (7 = 1 + {AC^) / {?,^''^'k). The stationary value of the 
BEG width is given by the fixed point of Eq. (jTTjl 

In the dimensional variables the BEG width is = [(27r)^/^/3^/®]V^/2. Linearizing near 
this solution the variational equation ()27|1 . we obtain the frequency of oscillations 2ujz that 
coincides with the hydrodynamic calculations based on the local field approximation (see 
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38, 



the next section). Equation ()27|) belongs to so-called Ermakov- Penney equations |l2 . 
and the solution is 

<t) = /^^ + ^' (28) 
where functions /? and c are linearly independent solutions to the equation 

yu + f{t)y = 0, (29) 

W = (3ct-c(3t is the constant Wronskian, and /5(0) = a(0), /?'(0) = a'(0), c(0) = 0, c'(0) = 1. 
For a periodic function /, Eq. is the so-called Hill equation which has been extensively 
studied j^]. So, in spite of the nonlinear character of Eq. the resonant response of gas 
has a linear character and, in particular, the frequency of oscillations does not depend on 
the amplitude. 

B. Hydrodynamic approach 

Eq. (j^Bjl can be cast in the form of Landau hydrodynamic equations by setting u = 



y/pexp[i6) ll|,|4l|. Furthermore we get the Thomas- Fermi solution by neglecting the kinetic 



energy term with respect to the interaction term, so that the equations read 

Pt = -{pv)x, (30) 
vt = -vvx - f(t)x + 2pp^, (31) 

where the velocity field v is defined by 9(t, x) = v{x', t)dx' . The equilibrium profile for p 
corresponds to a stationary solution of the form u{t, x) = exp{—ipt)u{x) where the chemical 
potential fi is related to the normalized number of atoms C = J \u\'^dx through the identity 
p = \piC In = N. The equilibrium profile has finite extension 



for X G (— -\/2a, \/2a), the velocity field is xb{t), and a and b satisfy the coupled equations 

at = ab, (33) 



bt = -b' - fit) + ^. (34) 
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Accordingly a satisfies the closed form equation 

att + f{t)a = (35) 

where C = jj? = 2(7^/7?^. Note that we have normalized the density profile so that a/\/2 is 
the rms width, which is the same as for the Gaussian ansatz We can thus compare the 
result (j27j) obtained with the variational approach using the Gaussian ansatz with the result 
(IH5|1 obtained with the hydrodynamic approach. Taking into account that the hydrodynamic 
approach is derived in the framework of a large number of atoms to neglect the kinetic term, 
we get that both approaches give the same effective equation, up to a small mismatch in the 
numerical value of C* ^ 0.245C^ (variational approach) and C ~ 0.203C^ (hydrodynamic 
approach). This departure originates from the fact that the two stationary profiles do not 
coincide. Note that the stationary point is C*^/^, so that the difference is around 5% which 
is negligible in practical situations. Eventually, the BEG dynamics is found to be governed 
by the same effective equation according to both approaches. 



C. Periodic modulations of the trap potential 

The dynamics of the Hill equation driven by a periodic modulation f{t) = 1 + /i sin(nt) 
is characterized by a parametric resonance phenomenon studied in j3^. In particular the 
stability of the solutions to the Hill equation depends on the parameters /i and VL. The 
theoretical prediction is that the stationary solution ag is unstable when the modulation 
frequency is close enough to the eigenfrequency of the BEG, more precisely 

|f]_2|<^. (36) 

We have performed numerical simulations of the PDE ()2fi|l to check this theoretical pre- 
diction. We have initiated the PDE with the initial condition given by the theoretically 
stationary Gaussian profile with C = 57r/\/2. The initial width is then Ug ~ 2.36. Note in 
Figure inithe presence of a cone of instability in the (/i, Q) landscape, in full agreement with 
Eq. ()36|). We have found numerically that the dynamics is unstable if |f2 — 2| < 0.55|/i|. 

Note that the equation for the oscillations of a cloud of nonineracting fermions in a 
time-dependent elongated trap has the same form as Eq. (f?fj) (see for example j^^). As a 
consequence the linear parametric resonance for the width of fermionic cloud at the same 
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FIG. 6: Contour levels of the oscillation amplitude of the BEC width in TG regime in presence of 
periodic modulation of the trap potential f{t) = 1 + /i sin(r2t). The gray area corresponds to the 
configurations where a blow-up of the solution has been numerically observed. 

frequency exists. In that sense we can say that the Fermi-Bose mapping still exists for the 
Tonks gas in an oscillating trap potential. 



D. Random modulations of the trap potential 

We examine in this section the effects of random modulations of the trap potential of 
the form /(t) = 1 + r]{t) where t] is a zero-mean stochastic process. The dynamics then 



exhibits stochastic resonance as shown for instance in j42l. |43| . Contrarily to the periodic 
case, we always observe an exponential growth of the oscillations of the BEC width, unless 
the random modulation has a vanishing power spectral density in the vicinity of the resonant 
eigenfrequency 2. More precisely, we get that a periodically oscillates between the values 
Gg ± ttoscit) with Gosc that grows exponentially as 

ao,c{t) ~ exp f ^ + ^W^tJ , (37) 

where = /o°° (?7(0)?7(s)) cos(2s)(is and Wt is a Brownian motion. If is a white noise 
{v{^)v{'t)) = 2a'^6{t), then we simply have ac = cr"^- The long-time behavior of a is dominated 
by the deterministic exponential growth term exp{act/4:) with very high probability because 
Wt ~ Vi. Note however that taking the expectation of Eq. (jHTj) yields a different exponential 
growth rate 

<«L>'^' it) - exp {^t) . (38) 
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This is due to the fact that some exceptional reahzations of the random fluctuations may 
induce very strong oscillations, and these exceptional realizations actually impose the value 
of the expected value. 

In the white noise case {v{0)r]{t)) = 2o'^6{t), using the same linearization procedure as 
in Section 1111 CI we can get precise expressions for the rms amplitude as long as stochastic 
resonance can be neglected 



, , , , aag / sin(2a;t) 
rms(a)(t) = — (39) 

where a; = 2 is the eigenfrequency. The simplest way to take into account stochastic reso- 
nance is to multiply the previous expression by the exponential damping term exp(o"^t/2). 
There exists a more accurate way based on Ito's calculus. Assume that the initial state is 
a(0) = ao, a'(0) = 0. The column vector X = ((a^) , {ab) , (&^))^ satisfies the closed system 

2 ^ 



^ = MX, M 
dt 



-1 1 
y2a2 -2 y 



(40) 



starting from X(0) = Vq = (^al,0,C /alj . As a first application we can compute the 
exact expression of the largest eigenvalue of M which governs the exponential growth of the 
modulation (a^): 

(54(t2 + 6^48 + 81(7^)^^^ - 12 



A. 



1/3 



(54(t2 + 6V48 + 81^4) 

Note that we recover formula ()38|) by expanding this expression for a ^ 1: Xmax — + 
O(cr^). A straightforward numerical integration of Eq. gives the exact evolution of (a^). 
We have performed numerical simulations of the PDE system (|2b|) with a random modulation 
of the trap potential to check the predictions obtained with the variational approach. We 
have taken the model where ri{t) is stepwise constant over elementary intervals with duration 
tc and takes random values uniformly distributed between —1 and 1. With tc = 0.06 we 
then have ac = 0.01. We compare the results of numerical simulations with the theoretical 
predictions in Figure We can see that (jHEl) efficiently predicts the initial growth of the 
oscillation amplitude, but it is necessary to take into account stochastic resonance when the 
amplitude becomes larger. 
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FIG. 7: Rms of the BEC width obtained from the averaging over a set of 1000 numerical simulations 
(solid line) and compared with the theoretical growth rate without stochastic resonance (dotted 
line) and with stochastic resonance (dashed line). The initial state is the Gaussian ansatz with 
C = 57r/\/2 and a(0) = flg = 2.36. 

VI. CONCLUSION 

In this work we have considered the resonances in collective oscillations of ID Bose gas 
under time-dependent variations of the trap potential and the effective nonlinearity. Two 
regimes have been studied - the mean field regime and the Tonks- Girardeau regime. The 
analysis shows that in the mean field regime the resonances are nonlinear and bistability 
exists in the vicinity of a critical frequency which is significantly above the eigenfrequency 
of the BEC. The dynamics is then characterized by stable oscillations with large amplitudes 
which depend on the frequency detuning between the frequency of the breathing mode and 
the modulation frequency. This type of dynamics is also predicted for the modulations of the 
trap as well as for the modulations of the atomic scattering length. In the Tonks- Girardeau 
regime the theory based on the nonlinear Schrodinger equation with quintic nonlinearity 
predicts the linear parametric resonance in the gas oscillations. The effect reflects the Bose- 
Fermi map existing for the Bose gas in this regime We also study the oscillations under 
random variations of the trap potential and effective nonlinearity. The analysis shows that 
the dynamics is nonlinear in the mean field regime. In the Tonks-Girardeau regime the study 
predicts stochastic parametric resonance. We also investigate the dynamics of ID Bose gas 
in the mean field regime under rapid and strong modulations of the atomic scattering length. 
This problem has recently attracted a great attention since the dynamically stable nonlin- 
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earity managed atomic matter solitons can be generated. We find the analytical expressions 
for the stationary value of the width and frequency of the slow secondary oscillations of 
the width. All theoretical predictions turn out to be well supported by direct numerical 
simulations of the ID GP equation and the quintic nonlinear Schrodinger equation. 
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